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Abstract. We use the modified quark-meson coupling and the quantum hadrodynamics models to study 
the properties of neutron star. Coupling constants of both models are adjusted to reproduce the same 
saturation properties. The onset of kaon condensation in neutron star matter is studied in detail over a 
wide range of kaon optical potential values. Once the kaon condensation takes place, the population of 
kaons increases very rapidly, and kaons become the dominant component, possibly making the neutron 
star matter a kaonic matter if the kaon optical potential is large. 

PACS. 97.60. Jd Neutron stars - 14.20.Jn Strange particles 



1 Introduction 

Observations of the masses of binary pulsars PQ , which are 
candidates of neutron stars indicate that the maximum 
mass of neutron stars are roughly around 1.5M . While 
the mass of a neutron star with only neutrons and protons 
is estimated to be about 2M Q , exotic degrees of freedom 
[2,3,4] other than the nucleons seem to reduce the maxi- 
mum mass of a neutron star close to the observed values, 
which implies that the exotic degrees of freedom such as 
creation of hyperons, Bose-Einstein condensation, strange 
matter, and quark deconfincment seem to be needed to 
reproduce the observed masses of neutron stars. 

In this work, we consider the strangeness degrees of 
freedom by including hyperon creation and kaon conden- 
sation in the neutron star matter. The equation of state 
(EoS) of dense nuclear matter is sensitive to the in-medium 
interaction of the hyperons and kaons. Recently, the mag- 
nitudes of the kaon-nucleus potential in matter have at- 
tracted much attention. Some calculations [SUHIE] show 
that the real part of the if "-nucleus optical potential Uk- 
is shallow (Uk- ~ ~ 50 MeV), but some other calculations 
suggest that U K - can be as large as about —120 MeV [HI 
E] or even close to —200 MeV [10]. We consider the pos- 
sibility of deep optical potential of kaons in nuclei and 
explore the consequences in the composition of neutron 
star matter and the mass-radius relation of the neutron 
star. 

To treat the dense nuclear matter we employ two dif- 
ferent relativistic mean field models; the modified quark- 
meson coupling (MQMC) model [11] and the quantum 
hadrodynamics (QHD) model [12]. In the MQMC model 
the nucleons and hyperons in the baryon octet are treated 
as MIT bags, whereas in QHD they are assumed to be 



point particles. The parameters of the two models are cali- 
brated to produce exactly the same saturation properties. 
By comparing the results from the two different models 
we can investigate the model dependence at high densi- 
ties. We find the onset densities of the kaon condensation 
and the compositions of matter at high densities are non- 
negligibly model dependent. The mass-radius relation of 
the neutron star can significantly differ depending on the 
interaction of the kaon in nuclear matter. 



2 Models 

The model Lagrangian consists of several terms for the 
octet baryons, exchange mesons, leptons and kaons; Ctot — 
£b + £m + £i+£k- In the mean-field approximation, each 
term reads 
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where B denotes the sum over the octet baryons and I 
stands for the sum over the free electrons and muons. The 
interactions between the non-strange light quarks (u and 
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d) are mediated by a, to and p mesons, and a* and <fi 
mesons are introduced to take into account the interac- 
tions between s quarks. One of the most important differ- 
ences between MQMC and QHD lies in the treatment of 
baryons; the former treats baryons as composite systems 
of quarks, and the latter assumes baryons as point parti- 
cles. As a result, the analytic form of the effective mass of 
a baryon differs greatly as follows. 

In MQMC, the effective mass of a baryon in matter 

m* B (a, a*) can be written as m* B = \J~E^ B ~- ~ EgT^fO^' 
where Eg is the bag energy of a baryon, R the bag radius 
and x q the eigenvalue of the quarks in the bag. The bag en- 
ergy of a baryon is given by E B = Y. q Ti^^r +^kR 3 B b , 
where Bg and Zb are the bag constant and a phenomeno- 
logical constant for the zero-point motion of a baryon 

B, respectively, and fl q — \j x \ + (R m q) 2 , where m*(= 

m q — g^a — g^,a*) is the effective mass of a quark. In the 
MQMC model, the bag constant Bb is assumed to depend 
on density, and we use the extended form of Ref . [3] . 

In QHD, the effective mass of a baryon is written as 
m* B = in b — gcrB<3 — gcr'BO* ■ We require MQMC and 
QHD to satisfy identically the same saturation proper- 
ties. To produce the compression modulus within a rea- 
sonable range one needs to include self-interaction terms 
of the er-field represented by [/Q HD = c 3 + \gs c 4 in 
the Lagrangian for QHD. The meson term in the QHD 
thus reads £^ HD = C M - t/? HD 

The kaon is treated as a point particle in both MQMC 
and QHD models in this work, and its effective mass is 
given as m* K = mK—goKO—g a »K&*- Covariant derivatives 
in Eq. (|U) include interactions with vector mesons through 
£> M = <9 M + ig^K^n - ig$K$» + i\g P KT ■ P^- A plane wave 
solution for the kaon field equation gives the dispersion 
relation of the anti-kaon 

U) K = m* K - g^K^a + 94,k4>o - 9 p k^Po3- (5) 

Nuclear saturation properties are relatively well known, 
and the values quoted in the literatures lie in the range 
0.15 ~ 0.17 fm -3 for the saturation density po, 15 ~ 16.3 
MeV for the binding energy per nucleon Eb, 30 ~ 35 MeV 
for the symmetry energy a sym , 200 ~ 300 MeV for the 
compression modulus K, and (0.7 ~ 0.8)mjv for the ef- 
fective mass of the nucleon. In this work we assume the 
saturation properties po = 0.17 fm~ 3 , Eb = 16.0 MeV, 
a sym = 32.5 MeV, K = 285 MeV, and m* N = 0.78m N . 
These conditions determine the coupling constants of non- 
strange mesons to u and d quarks in the MQMC, and to 
nuclcons in the QHD. These coupling constants are sum- 
marized in Table [1] 

For the coupling constants between the non-strange 
mesons (a, u>, and p) and hyperons we use the quark count- 
ing rule assuming that the s-quark does not couple to u 
and d quarks. This gives us g s a = = g s p = 0, and these 
relations determine the coupling constants of non-strange 
mesons to hyperons in both MQMC and QHD. For the 
coupling constants of strange meson to hyperons, we as- 



Table 1. Non-strange meson coupling constants for the 
MQMC (left four columns) and for the QHD (right four 
columns) models. g P N in QHD is the same as g q p in MQMC. 
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sume SU(6) symmetry, which gives us g J, = V2g^ d and 
9% = V2gZ' d - 

We have five meson-kaon coupling constants, ga-K, g^K, 
g P K, g<r*K and g<j>K- The quark counting rule is employed 
for g^K and g P K- g<r»K can be fixed from /o(980) decay, 
and g^x from the SU(6) relation V^g^K — g-n-np- g<j'K 
and g^K thus fixed are 2.65 and 4.27, respectively. The 
remaining coupling constant, g a Ki can be related to the 
real part of the optical potential of a kaon at the satu- 
ration density through Vr- = —{g a K<J + g^K^o)- The 
resulting g a x values for several C/jf- are given in Tabled 

Once the coupling constants arc determined, one can 
obtain the EoS of neutron star matter by solving self- 
consistently the equation of motion of exchange mesons, 
charge neutrality condition, baryon number conservation, 
and /3-equilibrium conditions for baryons and kaons. The 
mass-radius relation of the neutron star can be obtained 
by inserting the EoS into the Tolman-Oppenheimer-Volkoff 
(TOV) equation. 



3 Results 

A particle fraction is defined by the density of a parti- 
cle divided by the baryon density. Numerical results of 
the particle fraction is shown in Fig. [1] The onset den- 
sity of the kaon condensation p cr i t decreases as |Z7jf-| in- 
creases. This is expected because a larger g a x makes m* K 
smaller, which then causes the chemical equilibrium of 
kaons, lok = p n — P P fulfilled at lower densities, where 
A*n(p) is the chemical potential of the neutron (proton). 
An interestring feature is that once the kaon is created, 
the density of K~ increases very rapidly and dominates 
the particle population over the hyperons and even the 
nucleons. This can be partly attributed to the fact that 
the contribution from the w-meson to the octet baryons 
is repulsive whereas the contribution from the w-meson 
to the kaon is attractive. The u-meson term in the en- 
ergy of K~ in Eq. ([5|) has a negative sign and is thus 
attractive, but it is repulsive for the chemical potential of 
the octet baryons. The w-meson enhances the population 

Table 2. g aK for several U K - values in MQMC and QHD. 
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Fig. 1. The relative populations of particles in nuclear matter 
calculated from the MQMC and QHD models are displayed 
in the left and right columns, respectively, for U K - = —120, 
-140 and -160 MeV. 



of K~ but suppresses baryons, and thus the kaon den- 
sity increases rapidly. In addition, due to the competition 
between the negatively charged hyperons and K~ in the 
charge neutrality condition, the negatively charged hyper- 
ons are highly suppressed and in some cases are not even 
created at all as soon as the kaon condensation sets in. 
Positively charged hyperons, on the other hand, receive 
the opposite effects from the kaon condensation, and S + 
is created at lower densities as |£/if-| increases. The pro- 
ton density is also enhanced by large abundance of if - , 
which facilitates in turn the enhancement of E + popu- 
lation through the chemical equilibrium condition of the 
positively charged hyperons. 

We find that the onset density of kaon condensation 
pcrit from the MQMC model is lower than that from QHD. 
For U K - = -120, -140, and -160 MeV, p clit from MQMC 
are 5.9po, 3.8po and 3.0/9o, respectively, while they are 
9.8/9o, 4.3po an d 3.3po in QHD. The model dependence of 
Pcrit becomes less pronounced for a larger \Uk- I value. 

Kaon condensation leads to a phase transition from 
a pure hadronic phase to a kaon condensed one with a 
mixed state of hadron and kaon phases. To deal with such 
a mixed phase one needs to use Gibbs conditions. How- 
ever, for simplicity, we employ Maxwell condition for the 
mixed phase, which is often used as an approximation to 
the Gibbs condition. Though the Maxwell construction 
may not give us very accurate results particularly for the 
EoS, it can still provdie us with reasonable results for the 
mass-radius relation of the neutron star. 



Fig'- 121 shows the EoS including kaon condensation, but 
calculated without using either Maxwell or Gibbs con- 
dition. There is an instability region where pressure de- 
creases as the density increases. If such an instability exists 
in the interior of a neutron star, there is a high pressure re- 
gion on top of a low pressure one. The low pressure region 
cannot sustain the high pressure region, and the unsta- 
ble region will collapse, and the matter will be rearranged 
to fix the instability. Including such a rearrangement in 
the calculation of the EoS is equivalent to implmcnting 
Maxwell or Gibbs condition. With Maxwell condition, the 
unstable part is replaced by a flat line, which is deter- 
mined by the equal area method. (The flat lines are not 
plotted in Fig.[2]to avoid too many lines in one figure but 
should be obvious.) From Fig. [5]we can see that kaon con- 
densation begins at a point where the EoS curve deviates 
from the solid curve obtained without kaons. The kaon 
condensation makes the EoS softer than the one without 
it. As | £//£•- 1 becomes larger, the EoS becomes softer. By 
comparing the two solid curves in the left and in the right 
panels we see that the EoS without kaons is softer with 
the QHD model, but the effect of kaon condensation is 
more significant in the MQMC model. 

Once the EoS is determined, it is straightforward to 
solve the TOV equation to obtain the mass-radius rela- 
tion of a neutron star. Fig.[3]shows the mass-radius curves 
calculated by the MQMC (left) and QHD (right) mod- 
els. For Uk- = —120 MeV, the effect of kaon conden- 
sation is rather small and negligible particularly for the 
QHD case. The maximum mass of a neutron star with 
only nucleons is about 2M Q , as mentioned in the Intro- 
duction. When hyperons are created, the maximum mass 
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Fig. 2. The EoS without using any equilibrium condition in 
the mixed phase. The left figure is from MQMC, and the right 
figure is from QHD. The unstable region where the pressure 
decreases with increasing density has to be remedied by either 
Maxwell or Gibbs condition. 
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Fig. 3. The mass-radius relation from the MQMC (left) and 
QHD (right) models. Both figures are obtained with Maxwell 
construction. 
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Fig. 4. The particle fraction with repulsive in-medium inter- 
action of E hyperon in the QHD model. Us is fixed as +30 
MeV, and U K - is chosen as -120 (left) and -160 (right) MeV. 



of a neutron star decreases about 0.5M©. As the depth 
of optical potential for K~ increases, the effect of kaon 
condensation becomes more significant. If U K - = —160 
MeV, the maximum mass decreases by about 10 ~ 12% 
from the maximum mass without kaon condensation. In 
the case of QHD the maximum mass with Uk- = —160 
MeV becomes so small as 1.32M Q that most of the ob- 
served masses in Ref. [1] exceeds this (maximum) value. 
This raises a question whether a very deep optical poten- 
tial of K~ can be comptible with the existing observations 
of the neutron star mass. We will discuss this again in the 
next section. 

One of the unsettled issues in hypernuclear physics 
is the interaction of S hyperon in nuclear medium. Use 
of plain quark counting rule for the hyperon-meson cou- 
pling constants gives us attractive optical potentials of 
S hyperons in the range —40 ~ —30 MeV at the satu- 
ration density. Some studies of Z'-hypernuclei, however, 
indicates that the interaction of 27 hyperon in medium is 
repulsive and its optical potential at the saturation den- 
sity takes a positive value [13 J. To take into account this 
possibility, we have repeated the calculations using the 
QHD model with the repulsive 27-hyperon optical poten- 
tial of +30 MeV at the saturation density and show the 
results in Fig. 01 The particle compositions are displayed 
for Uk- — —120 and —160 MeV. By comparing these two 
figures with the figures on the right panels of Fig. [TJ one 
can see that the repulsive 27-hyperon interaction makes 
big changes in the population of hyperons, but the onset 
densities and the population of kaons are little affected. 
The maximum masses of a neutron star with Us = +30 
MeV are 1.52M©, 1.48M and 1.32M for U K - = -120, 
— 140 and —160 MeV, respectively. These results deviate 
from those with attractive 27-hyperon interaction only by 
about 1.4% at most. 



4 Discussion 

Masses of some pulsars fT4l[T5] that are more recently mea- 
sured are reported to lie much above the canonical value 
of Ref. pp. Some of the large masses are in the range of 
(1.7 — 2.1) M®, though the results depend on the binary 
systems. These values are even larger than the maximum 
mass of a neutron star composed of only nucleons and hy- 
perons (in short, a baryon star). If kaon condensate states 



are added to the baryonic state, it will soften the EoS and 
reduce the maximum mass further down. Then one con- 
fronts the question: Are exotic states such as kaon conden- 
sation ruled out by large mass neutron stars? To answer 
the question, we have to understand better the interac- 
tion of particles and state of matter at densities above 
the nuclear saturation. For instance, we don't know yet 
to what densities our mean field lagrangian models are 
valid. If the large mass neutron star does exist, we have to 
have a mechanism which makes the EoS at high densities 
stiffer than now. Introducing hard core repulsion between 
nucleons can make the EoS stiff at high densities. The 
effect of hard core has been considered in the high den- 
sity/temperature limit region [T6lfT7lfT8] . and it is shown 
that it makes the deconfined quark phase the most sta- 
ble state of matter at high density and/or temperature. 
One can take into account such effects in the neutron star 
matter together with exotic states. The investigation is in 
progress. 
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